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Bamboo, as a natural hierarchical cellular material, exhibits remarkable mechanical properties including 
excellent flexibility and fracture toughness. As far as bamboo as a functionally graded bio-composite is 
concerned, the interactions of different constituents (bamboo fibers; parenchyma cells; and vessels.) 
alongside their corresponding interfacial areas with a developed crack should he of high significance. Here, 
by using multi-scale mechanical characterizations coupled with advanced environmental electron 
microscopy (ESEM), we unambiguously show that fibers' interfacial areas along with parenchyma cells' 
boundaries were preferred routes for crack growth in both radial and longitudinal directions. Irrespective of 
the honeycomb structure of fibers along with cellular configuration of parenchyma ground, the hollow 
vessels within bamboo culm affected the crack propagation too, by crack deflection or crack-tip energy 
dissipation. It is expected that the tortuous crack propagation mode exhibited in the present study could be 
applicable to other cellular natural materials as well. 

As one of the most renewable resources on Earth, bamboo is the fastest-growing and highest-yielding 
natural cellular material, which reaches maturity within months and ultimate mechanical properties 
within few years 1 . Well known as a functionally graded hierarchical bio-composite, bamboo comprises 
three fundamental tissues named epidermis, vascular bundles and parenchyma ground. The thick epidermis is the 
shell of the bamboo whereas the vascular bundles are the longitudinal tissues supporting the whole bamboo, with 
the ground parenchyma occupies the rest of the organ. Within each vascular bundle, the role of vessels and 
phloem is to transport water and nutrients 2 whereas they are all surrounded by fibers 3 (see Fig. 1). In terms of 
volume fraction, the fibers and cellular parenchyma form the majority of bamboo culm (for Phyllostachys edulis 
species: ~40-60% 4,5 and ~20-60% 6,7 , respectively, depending on location, local climate, age...); whereas vessels 
and phloem make up the reminder. In view of the weight-to-weight basis, tensile strength, Young's modulus, 
compressive strength and interlaminar shear of bamboo is reasonably comparable with conventional structural 
materials such as low carbon steel and fiber glass reinforced plastics 8 . So far, the interesting mechanical properties 
of bamboo have been mainly attributed to the presence of fibers within the bamboo culm 81 ". However, for such a 
complicated hierarchical structure, one wishes to understand the role of other structural features, such as 
parenchyma cells and hollow vessels, on the mechanical performance of bamboo. And getting a deep under- 
standing on bamboo's hierarchical features, particularly at cellular level, can be quite helpful in designing 
biomimetic polymeric, metallic composites. 

Although a lot of attempts have been made, so far, to investigate the functionally graded structure along with 
bulk properties of bamboo 4,5,814 , very few attempts have been made to investigate the crack growth mode along 
with fracture mechanisms at micro-scale or cellular level. Among those earlier efforts, Shao et. al. 14 explored the 
behavior of interlaminar fracture using double cantilever beam specimens, and basically illustrated that the crack 
propagation develops along the longitudinal interface between the fibers and ground tissue. In another attempt by 
Low et. al. 8 , the excellent damage tolerance of bamboo was attributed to the interlay and simultaneous presence of 
crack deflection, fiber debonding and crack bridging as the major energy dissipative processes. Similarly, Tan et. 
al. 10 demonstrated that, in the course of bending deformation on a single edge notched specimen, the crack growth 
occurs by deflection into interlaminar boundaries. Pertaining to the studies conducted earlier, it is likely that the 
interaction between a developed crack with functionally graded fibers has been paid great attention to; whereas 
the role of cellular parenchyma ground along with presence of hollow vessels within the bamboo culm have been 
largely neglected. However, in light of considerable volume fraction of parenchyma cells along with hollow vessels 
(—57.3 ± 2.5% and —5.6 ± 0.8%, respectively, from our samples' microstructure analysis, see Methods), it is 
reasonable to speculate that the bamboo's remarkable mechanical behavior could be also stemmed out from the 
cellular configuration of parenchyma ground along with the possible crack interaction with hollow vessels. So in 
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the present study, we will focus on investigating the interactions of 
not only bamboo's fibers but also parenchyma ground alongside 
hollow vessels with developed cracks, in a holistic approach, to obtain 
a comprehensive understanding of the underlying mechanisms for 
bamboo's superior fracture toughness 15 . 

Results 

Microhardness indentation-induced crack growth. To investigate 
the interaction of different constituents of bamboo (see Fig. 1) with a 
developed crack, a crack was intentionally initiated by microhard- 
ness indentation along the radial direction (RD) of bamboo culm 
(see Fig. 2). Our microstructural investigations conducted on an 
indentation-induced crack (see Fig. 3a) clearly demonstrated that, 
once the initiated crack reached the critical condition for further 
crack growth, it propagated along the parenchyma cells' bound- 
aries (see Figs. 3b-c and i). The exhibited crack growth mode 
continued as long as the crack was approaching the fibers' bundles. 
Similarly, within individual fibers' bundles, the crack propagated 
further through the interfacial areas (see Figs. 3d-e and i). 
Pertaining to the exhibited microhardness-induced crack growth 
along the RD (Fig. 3), it was shown that the crack was mostly 
deflected by an average deflection angle (0) of 9 RD ~ 25°-60° and 
Bro ~ 45°-60° owing to the occurrence of interfacial fractures within 
the cellular parenchyma matrix and honeycomb fibers' bundles, 
respectively. Investigations further revealed that, in addition to 
cellular parenchyma matrix along with honeycomb fibers' bundles, 
the hollow vessels within bamboo culm could additionally affect the 
crack growth (see Fig. 3a): as cases demonstrated in Figs. 3f-h, the 
hollow vessel clearly deflected the crack and dissipated its driving 
energy. In fewer cases where cracks were not deflected and eventually 
entered into the hollow vessels, the vessels were essentially func- 
tioning as a crack-tip energy absorber. In either scenario (Fig. 3i), 
the vessels could effectively enhance the fracture toughness, not 
mentioning that the considerable volume fraction of hollow vessels 
also reduce the overall weight of bamboo. 

Tension-induced crack growth. To further validate the exhibited 
crack growth mode in the longitudinal direction (LD) of bamboo 
culm, tensile tests were conducted on dog-bone shaped samples 
(without pre-notch), prepared from bamboo's longitudinal sec- 
tions (see Fig. 2a) to assess the interaction of different constituents 
with a tension-induced crack. Microstructural analysis conducted on 
the tensile deformed samples, primarily, revealed that fibers' pull out 
(see Fig. 4a) was the most prominent feature in tensile fracture 
surfaces which could be simply correlated to the exhibited crack 
growth along the fibers' interfacial areas in both radial and 
longitudinal directions. More importantly, as displayed in Fig. 4b, 
the presence of those intact parenchyma cells was also regarded as a 
dominant feature in tensile fracture surfaces, which could be 
naturally attributed to the interfacial fracture within the cellular 
parenchyma ground, in the course of tensile loading. Similar to exhi- 
bited crack deflection along RD during microhardness indentation, 
the crack deflection angle (Old) along the LD, owing to tension- 
induced interfacial fracture, was estimated to be almost —90° in 
average (see Fig. 4c) within both parenchyma matrix and fibers' 
bundles here. 

In view of the compatibility amongst the interaction of different 
constituents of bamboo with both indentation/tension-induced 
cracks (see Figs. 3-4), it can be concluded that the cellular par- 
enchyma ground along with honeycomb fibers accommodate the 
crack growth via their interfacial boundaries, and hollow vessels also 
contribute in changing the crack growth mode via either absorbing or 
deflecting the crack. Likewise, pertaining to the observed crack 
growth modes in the case of both indentation and tensile loading 
(see Fig. 3i and Fig. 4c), the crack deflection and propagation in real 
three-dimensional (3D) space can be therefore speculated: the 



indentation-induced crack exhibits the tortuous crack growth in 
bamboo's radial direction (RD); whereas the tension- induced crack 
demonstrates the crack growth in both radial (RD) and longitudinal 
(LD) directions. Together we have obtained a complete picture 
showing the interactions of bamboo's different constituents with a 
developed crack in a real bamboo culm structure. 

Discussion 

It is commonly known that, as a functionally graded bio-composite, 
the toughening of bamboo mostly happens during the crack growth 
and not in the course of crack initiation 8,913 . Hence, the more difficult 
is the crack growth, the higher is the fracture toughness. To date, the 
exhibition of bamboo's remarkable mechanical properties has 
mainly been attributed to the presence of fibers within the bamboo's 
culm 8,9,13 . To get a quantitative understanding of the property differ- 
ence among bamboo's different phases, we also performed micro- 
and nano-mechanical characterizations on bamboo's individual 
constituents and confirmed that, fibers are indeed the strongest 
phase (elastic modulus Ef = 22.8 ± 2.8 GPa by nanoindentation 
and E f = 30.1 ± 3.0 GPa by micro-tensile testing; tensile strength 
c>f = —1000 ± 300 MPa; see Methods) whereas parenchyma cells 
are comparatively weaker (elastic modulus E p = 3.7 ± 0.4 GPa by 
nanoindentation). However, our experimental observations on frac- 
tured samples strongly suggested that the cellular configuration of 
parenchyma ground along with presence of hollow vessels within the 
bamboo culm (with their considerable volume fractions ~ 57.3 ± 
2.5% and —5.6 ± 0.8%, respectively) should be also responsible 
for the bamboo's remarkable fracture toughness (averagely — 
56.8 MPa.nv) 15 : the exhibited interfacial fracture within the cellular 
parenchyma matrix along with honeycomb structure of fibers 
revealed that the interfaces of parenchyma cells as well as that of 
fibers were sufficiently weaker compared to the fracture toughness 
of both parenchyma cells and fibers themselves, respectively. 
Pertaining to the demonstrated interfacial fracture, the high fracture 
toughness of bamboo could be partly justified. The prevailed tor- 
tuous crack growth mode effectively increased the final fracture sur- 
face which consequently amplified the driving energy required. 
Likewise, the remarkable fracture toughness could be also attributed 
to the presence of hollow vessels within the bamboo culm as they 
could either deflect the crack growth direction or absorb the crack-tip 
energy (see Fig. 3a and 3i). Hollow vessels should be of high signifi- 
cance as their presence greatly reduces the overall weight of bamboo, 
whereas in view of the geometry and configuration, their function as 
a stress raiser in the course of loading is unlikely. 

In view of the exhibited tortuous crack growth mode both in RD 
(see Fig. 3i) and LD (see Fig. 4c), the reduction in driving force for 
crack propagation, owing to both occurrence of interfacial fracture 
and crack tips' interaction with hollow vessels, can be likewise 
justified by a geometry- dependent model proposed by Hanlon 
and Suresh et. al. 1617 for polycrystalline materials: considering the 
exhibited crack growth mode, it is assumed that the cellular bam- 
boo structure was mostly experiencing the mode I crack opening in 
the course of loading 14 and any deviation of the crack away from 
the mode I growth plane (straight crack growth) led to a mixed 
loading mode at the tip of the crack 16 . In the case of the exhibited 
tortuous crack growth in both radial (RD) and longitudinal (LD) 
directions, it is assumed that the crack kinked at an average angle 
inclined 8 with respect to the mode I plane. So, the local effective 
stress intensity factor at the tip of the crack could be then approxi- 
mated by 16 : 

Keff=COS 2 QKl (1) 

According to the equation, comparing to the case of having a 
straight crack growth (Ki), the crack deflections, owing to either 
occurrence of interfacial fracture among parenchyma cells (6 RD — 
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Figure 1 | SEM micrographs of the raw bamboo culm with different constituents, zoom-in views of bamboo's vascular bundles along with the 
parenchyma ground and bamboo fibers along the transversal ((a), (c), (e) and (g)) and longitudinal ((b), (d), (f) and (h)) directions. As displayed, fibers 
and parenchyma cells, comparably, possess the majority of bamboo culm whereas vessels possess less contribution. 



25°-60°; 9 LD ~ 90°) and fibers (9rd ~ 45°-60°; 9 LD ~ 90°) or further crack propagation 16 . This can be simply realized as the 
interaction with hollow vessels (6 varies according to the vessel effective stress intensity factors (K eff ) owing to crack deflections 
curvatures), would apparently reduce the overall driving force for along the RD within the corresponding constituents (K c ff (RD)), 
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Tensile Loading 



(a) 



(b) 




Figure 2 | Samples and experimental configuration, (a) microhardness (left, cubic) and tensile (right, dog-bone shaped) testing samples, along with the 
microhardness indentation configuration; (b) schematic representation of microhardness probe loading direction (the red arrow) with respect to a 
vascular bundle. 



which were estimated to be K^ ff (RD) ~ (0.75 — 0.95)K| and 
K^ ff (RD) ~ (0.75-0.85)K[, for the parenchyma matrix (K p ) alone 
and fiber bundles (K r ) alone, respectively. Similarly, along the LD, 
the effective stress intensity factors were estimated to be 
K^ ff (LD) ~ (0.5)Kj and K^ ff (LD) ~ (0.5)K| for the parenchyma 
matrix (K p ) and fiber bundles (K*), respectively. The reduction in 
overall effective stress intensity factor owing to the crack deflection 
by vessels could be likewise justified in a similar manner. Despite 
the different stress intensity factors (K e ff) of bamboo's respective 
phases owing to structural heterogeneity, unlike polycrystalline 
solids, the combined effect according to aforementioned estima- 
tions and their interplays will result in an overall reduction of the 
driving force for crack propagation in the whole bamboo structure. 
Therefore, as demonstrated experimentally and validated numer- 
ically, the occurrence of interfacial fractures within the cellular 
parenchyma matrix as well as honeycomb fibers' bundles could 
be considerably responsible for bamboo's remarkable fracture 
toughness, along with the contribution from the hollow vessels. 
Lastly, it might be noted that, such interactions amongst the struc- 
tural features and a propagating crack has been likewise demon- 
strated for other plant tissues: for example, in Yew and Spruce 
wood, fibers with different shapes and sizes retard the crack growth 
via fiber bridging whereas the holes within the structure arrest the 
crack in the course of its propagation 18 , while in the case of Green 
Sapwood (Pinus Sylvestris), crack propagation often occurs 
through the middles lamella of earlywood, causing deflection 
around the cells, followed by arresting as it approaches the late- 
wood layers 1 ". 



Conclusions. In summary, the remarkable mechanical behavior of 
bamboo is stemmed out from the hierarchical configuration of 
parenchyma ground as the matrix along with fibers as the 
reinforcement. The occurrence of interfacial fractures within the 
cellular parenchyma matrix along with honeycomb structure of 
fibers in both radial and longitudinal directions is mainly 
responsible for the remarkable fracture toughness of bamboo by 
reducing the overall effective stress intensity factor. Likewise, the 
hollow vessels within the bamboo culm can further tailor the 
fracture properties by deflecting crack growth or absorbing crack- 
tip energy. Therefore a complete picture showing the tortuous crack 
propagation mode in bamboo's hierarchical cellular structure has 
been unambiguously disclosed. It is expected that the present 
experimental findings concerning the bamboo's unique structural 
features versus its mechanical behavior shall be quite helpful in 
creating biomimicking structural designs along with the develop- 
ment of futuristic advanced composites. 



Methods 

Sample preparation and microstructure characterization. The moso or mao zhu 
bamboo (Phyllostachys edulis species), collected from raw bamboo plantations 
located in Jiangsu and Zhejiang provinces in China, were selected for this study. The 
samples were all mature bamboo (—5 years old) freshly cut from the middle section of 
stalk, and were kept at room temperature of ~23°C and relative humidity of —55- 
65% (reasonably close to the humidity level of moso bamboo's natural habitat) to 
avoid drying artifacts. Microstructural characterizations were conducted on the 
polished bamboo samples taken from both radial/transverse and longitudinal 
sections, to reveal the microstructure of different constituents. The characterizations 
were further continued to investigate the crack growth modes in the case of 
microhardness indentation-induced cracks along with fracture surface analysis after 
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Figure 3 | Microhardness indentation-induced crack growth in bamboo culm (a), along with: (b-h) high magnification micrographs representing the 
crack propagation within the corresponding constituents ((b-c): parenchyma cells; (c-e): fibers; (f-g): a hollow vessel; (g-h): fibers); (i) schematic 
representation of the indention-induced crack growth within different constituents of bamboo in radial direction (RD). 
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Figure 4 | Crack propagation in bamboo culm's three-dimensional (3D) 
structure revealed by a tensile fractured sample: (a-b) ESEM micrographs 
displaying the tensile fracture surfaces of the bamboo culm at high 
magnification, with focus on the (a) fibers pull-out and the (b) intact 
parenchyma cells; (c) schematic representation of crack growth within 
bamboo's different constituents from the view of longitudinal direction 
(LD). 

tensile deformations. For this purpose, a Philips™ XL30 FEG Environmental 
Scanning Electron Microscope (ESEM) along with a FEI Quanta IM 450 FEG Scanning 
Electron Microscope working at ESEM mode were used. The captured micrographs 
were also processed by TalyMap™ Universal image analyzer software to quantify the 
volume fractions of bamboo's different constituents. 

Mechanical characterizations. To investigate the interactions of different 
constituents alongside their corresponding interfacial areas with a developed crack, 
cracks were intentionally initiated along the radial direction (RD) of bamboo culms 
by a Fischerscope 1M HM2000 microhardness tester. For this purpose, a small cubic 
sample was subjected to microhardness indenter with an indentation load of 
2000 mN (see Fig. 2a-b). Later on, to further validate the observed indentation- 
induced crack growth mode, tensile experiment was conducted on an intact dog-bone 



shaped sample, prepared from bamboo's longitudinal direction (LD, see Fig. 2a), by 
using a MTS Alliance™ RT/30 Material Testing System. Similar to indentation- 
induced crack, pulling the dog-bone bamboo samples along the LD opened up a 
fracture surface along the RD which subsequently led to the crack propagation and 
fracture failure. Lastly, to obtain a quantitative understanding of bamboo's different 
constituents' mechanical properties, nanoindentation and micro-tensile tests were 
also performed. For nanoindentation experiments, a Hysitron TI 750 Ubi™ 
nanoindenter with a calibrated cube-corner tip (curvature radius — 100 nm) was 
employed. Load control experiments were used to indent the polished surface of 
individual fibers along with parenchyma cells, respectively. For micro-tensile tests on 
individual fibers, a Gatan Microtest 1M 200 tensile tester with 2 N load cell (resolution 
0.001N) was used under displacement control mode at a constant rate of 0.1 mm/ 
min. All mechanical characterizations were conducted on as received samples at 
ambient laboratory environment (temperature 23°C and relative humidity ~ 50- 
60%, respectively), to simulate natural bamboo's realistic deformation conditions. To 
further validate the crack propagation mode in natural bamboo with even higher 
moisture content, additional experiments for bamboo with high humidity (~75- 
85%) were performed (results shown in the Supplementary Information), by micro- 
hardness indentation experiments on pre -moisturized samples from high-humidity 
environment (relative humidity 90%, kept for overnight, >12 hours). 
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